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Abstract 

A new class of exact solutions of Einstein's field equations with a 
perfect fluid source, variable gravitational coupling G and cosmological 
term A for FRW spacetime is obtained by considering variable deceleration 
parameter models for the universe. The nature of the variables G(t), A(t) 
and the energy density p{t) have been examined for three cases : (i) 
exponential (ii) polynomial and (iii) sinusoidal form. The special types 
of models for dust, Zel'dovich and radiating universe are also mentioned 
in all these cases. The behaviour of these models of the universe are also 
discussed in the light of recent supernovae la observations. 
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1 Introduction 

In Einstein's General Relativity (GR), there are two parameters, the gravita- 
tional constant G and the cosmological constant A which are considered to be 
time- independent jl] . A possible time variable G was suggested by Dirac |2] and 
this has been extensively discussed in literature Since its introduction, its 

significance has been studied from time to time by various workers |S] — [7] ■ In 
modern cosmological theories, the cosmological constant remains a focal point 
of interest. A wide range of observations now compcUingly suggest that the uni- 
verse possesses a non-zero cosmological constant 8 . In the context of quantum 
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field theory, a cosmological term corresponds to the energy density of vacuum. 
The birth of the universe has been attributed to an excited vacuum fluctua- 
tion triggering off an inflationary expansion foUowed by the super-cooHng. The 
release of locked up vacuum energy results in subsequent reheating. The cosmo- 
logical term, which is measure of the energy of empty space, provides a repulsive 
force opposing the gravitational pull between the galaxies. If the cosmological 
term exits, the energy it represents counts as mass because mass and energy are 
equivalent. If the cosmological term is large enough, its energy plus the matter 
in the universe could lead to inflation. Unlike standard inflation, a universe with 
a cosmological term would expand faster with time because of the push from 
the cosmological term 0. Some of the recent discussions on the cosmological 
constant "problem" and on cosmology with a time-varying cosmological con- 
stant by Ratra and Peebles jJUj, Dolgov and Sahni and Starobinsky 
jl4| point out that in the absence of any interaction with matter or radiation, 
the cosmological constant remains a "constant". However, in the presence of 
interactions with matter or radiation, a solution of Einstein equations and the 
assumed equation of covariant conservation of stress-energy with a time- varying 
A can be found. This entails that energy has to be conserved by a decrease 
in the energy density of the vacuum component followed by a corresponding 
increase in the energy density of matter or radiation [see also Weinberg ^] and 
Carroll, Press and Turner [T^]. 

Recent observations by Perlmutter et al. ^7] and Riess et al. il8| strongly 
favour a significant and a positive value of A with magnitude K{Gh/c^) « 
2^0-123^ Their study is based on more than 50 type la supernovae with red- 
shifts in the range 0.10 < z < 0.83 and these suggest Friedmann models with 
negative pressure matter such as a cosmological constant (A), domain walls or 
cosmic strings (Vilenkin Garnavich et al. HOI)- Recently, Carmeli and 
Kuzmenko j21j have shown that the cosmological relativistic theory predicts 
the value for cosmological constant A = 1.934 x 10^'^^s^^. This value of "A" 
is in excellent agreement with the recent estimates of the High-Z Supernova 
Team and Supernova Cosmological Project (Garnavich et al. |2D|; Perlmutter 
et al. ^7]; Riess et al. JSI; Schmidt et al. These observations suggest on 

accelerating expansion of the universe. 

Since Dirac first considered the possibility of a variable G, there have 
been numerous modifications of general relativity to make G vary with time 
|23|. Recently a modification linking the variation of G with that of A has been 
considered within the framework of general relativity by a number of workers 
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The Einstein's field equations are a coupled system of non-linear differential 
equations. We seek physical solutions to the field equations for their applica- 
tions of relevance to cosmology and astrophysics. In order to solve the field 
equations we normally assume a form for the matter content and killing vector 
symmetries by space time Solutions to the field equations may also be 
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generated by applying a law of variation for Hubble's parameter proposed by 
Berman [SJ. In simple cases the Hubble law yields a constant value of decel- 
eration parameter. It is worth observing that most of the well-known models 
of Einstein's theory and Brans-Deke theory with curvature parameter k = 0, 
including inflationary models have a constant deceleration parameter as borne 
out by the studies of several authors [HI — ISHl- But redshift magnitude test has 
had a chequered history. During the 1960s and the 1970s, it was used to draw 
very categorical conclusions. The deceleration parameter was then claimed to 
lie between and 1 leading to the suggestion that the universe is decelerating. 
But today, the view is entirely different. Observations (Knop et al. |37]) Riess 
et al., of Type la Supernovae (SNe) allow us to probe the expansion history 
of the universe leading to the conclusion that the expansion of the universe is 
accelerating. So we can consider cosmological models with variable cosmologi- 
cal term and deceleration parameter. Vishwakarma and Narlikar |39) and Virey 
et al. 00] have reviewed the determination of the deceleration parameter from 
Supernovae data. 

Recently Pradhan and Otarod 41 have studied a model of the universe 
with time dependent deceleration parameter and A-term in presence of perfect 
fluid. This study has prompted us to focus upon the problem of establishing 
a formalism for understanding the relativistic evolution for a time dependent 
deceleration parameter in an expanding universe with variable G and A in pres- 
ence of perfect fluid. This paper is organized as follows. We present in Section 2 
the metric and the fleld equations. In Section 3 we deal with a general solution. 
In Sections 4, 5 and 6 we deal with the solutions in exponential, polynomial and 
sinusoidal form. Finally in Section 7 we mention our main conclusions. 



2 The metric and field equations 



In standard coordinates (a;°) = {t,r,9,(p) the Robertson- Walker line element 
has the form 



dr' 



1 — kr'^ 



+ r2(d6|2 + sin2 



(1) 



where S{t) is the cosmic scale factor. Without loss of generality the constant k 
takes only three values: 0, 1 or —1. The constant k is related to spatial geometry 
of a 3-dimensional manifold generated by t = constant. The Robertson- Walker 
spacetimes are the standard cosmological models supporting observations. For 
the case of variable cosmological term K{t) and gravitational constant G{t) the 
Einstein field equation 

Gab + Agab = 87TGit)Tab (2) 



yields 
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+ k) = SttGp + A, 



(3) 
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2^ + i^l±^^-8nGp + A, (4) 

for the line element ^ where a dot indicate partial differentiation with respect 
to t. From Eqs. (j3Jl and I^J we obtain the generalized continuity equation 

P + 3f(p + p) + §M+^ = 0. (5) 

This reduces to the conventional continuity equation when A and G are con- 
stants. In an attempt to obtain solutions to the field equations we assume, as is 
often done, that the classical conservation law, T.'^''b — 0, also holds. Then Eq. 
© leads to 

P + 3|(p + p), (6) 

SirpG + A = 0, (7) 

which facilitate the solution of the field equations. The result ^ is just the 
conventional continuity equation, and Q simply relates G and A and does not 
explicitly contain the scale factor S{t). 

3 Solution of the Field Equations 

Equations © - O are two independent equations in five unknowns S, p, p, G, 
and A. In this paper we first assume an equation of state 

P = 7P,0<7<1. (8) 

Secondly, on the basis of supernovae searches, we consider the deceleration pa- 
rameter to be variable and set 



SS H + H^ 



^(variable), (9) 



where H = S/S is the Hubble parameter. The above equation may be rewritten 
as 

--f6^=0. (10) 
The general solution of Eq. H1U|) is given by 

el i-^^dS = t + m, (11) 

where m is an integrating constant. 

In order to solve the problem completely, we have to choose J ^dS in such 
a manner that Eq. be integrable. 
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Without loss of generality, we consider 

J = In L{S), (12) 
and obtain from Eqs. ((TT|l and ((T^ 

L{S)dS^t + m. (13) 



Of course the choice of L{S), in Eq. (|13|l . is quite arbitrary but, since we are 
looking for physically viable models of the universe consistent with observations, 
we consider the following three cases: 

4 Solution in the Exponential Form 

Let us consider L{S) — -j^, where fci is arbitrary constant. 

In this case on integration of Eq. 113(1 gives the exact solution 

S{t) = fc2e"^*, (14) 

where ^2 is an arbitrary constant. 

Using Eq. Q in the differentiated form of Q, we obtain 

8.Gp=-^. (15) 
From Eqs. ©, Q and ©, we have 

2k 

8^(l+7)Gp=— . (16) 
From Eqs. (|15|l and 1(16(1 . we obtain 

^ = -3fci(l + 7), (17) 

which on integration gives 

p = ^e-3fci(i+7)t^ (18) 

where A is an integrating constant. Using Eqs. ((14(1 and ((18(1 in l(16() we get the 
value of G 

G^-r-r. rTT2- (19) 

47r(l + 7)Afc| ^ ' 



From Eq. J^Jl or (Q, we find A as 

"(1 + 7)fc2^ 



^ . " + 3fc^ (20) 
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From Eq. ifT^ . since scale factor can not be negative, we find S{t) is positive if 
k2 > 0. So in the early stages of the universe, i.e., near t = 0, the scale factor 
of the universe had been approximately constant and had increased very slowly. 
Some time later, the universe had exploded suddenly and expanded to a large 
scale. This picture is consistent with Big Bang scenario. 

From Eq. lfTH|) . we note that p{t) is a decreasing function of time and p > 
for all times. From Eq. H19|l. we observe that G is an increasing function of 
time. When the universe is required to have expanded from a finite minimum 
volume, the critical density assumption and conservation of energy-momentum 
tensor dictate that G increases in a perpetually expanding universe In 
most variable G cosmologies 03 B2| G is a decreasing function of time. But 
the possibility of an increasing G has also been suggested by several authors |44| . 

The behaviour of the universe in this model will be determined by the cos- 
mological term A ; this term has the same effect as a uniform mass density 
Peff = — A/47rG', which is constant in space and time. A positive value of A 
corresponds to a negative effective mass density (repulsion). Hence, we expect 
that in the universe with a positive value of A, the expansion will tend to accel- 
erate; whereas in the universe with negative value of A, the expansion will slow 
down, stop and reverse. From Eq. [20], we see that the cosmological term A is a 
decreasing function of time and it approaches a small positive value as time in- 
crease more and more. Recent cosmological observations (Garnavich et al. |19|. 
Perlmutter et al. Riess et al. Schmidt et al. [21]) suggest the existence 
of a positive cosmological constant A with the magnitude K{Gh/c^ k, 10"^^'^. 
These observations on magnitude and red-shift of type la supernova suggest 
that our universe may be an accelerating one with induced cosmological density 
through the cosmological A-term. Thus our model is consistent with the results 
of recent observations. In the following, we consider the variation of p{t), K{t) 
and G{t) with time for three different special cases. 

The expressions for p{t), K{t) and G{t) corresponding to 7 = 0,1,1/3 are 
given in Table 1. 
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A(t) 


Git) 









3fc2 


fee"!' 






1 










1 
3 




2ki ^ 




IG-r^Ak'^ 



Table 1: Values of pit), A(t) and Git) for dust, Zel'dovice and radiation 
exponential models. 

From Table 1, we observe that in all three models p > for all times and 
is a decreasing function of time. The cosmological term A is also decreasing 
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function of time and it approaches to a small positive value with large value of 
time. Thus our models have good agreement with recent observations. 



5 Solution in the Polynomial Form 



Let L{S) 



1 



, where /ca and are constants. 



In this case, on integrating, Eq. H13(l gives the exact solution 



S{t) — ait^ + + az, 



(21) 



where ai, a2 and are arbitrary constants. 

By using the same techniques as described in previous section 4, we find the 
expressions for p(t), G{t) and A(t): 



where B is an integrating constant. From Eq. H21f) . since scale factor can not 
be negative, we find S{t) is positive if ai, a2 and as are all greater the zero. 
From Figure 1, it can be seen that in the early stages of the universe, i.e., near 
t = 0, the scale factor of the universe had been approximately constant and had 
increases very slowly. At a specific time the universe has exploded suddenly and 
it has expanded to large scale. This fits with big bang scenario. 

From Eq. |j22Jl, we observe that p{t) > for i? > 0. We also see that p{t) is 
a decreasing function of time. Figure 2 clearly shows this behaviour of p{t). 

From Eq. H23|l . we note that G{t) is an increasing function of time. Figure 
3 also shows this behaviour of G as increasing function of t. From Eq. (|26|l . 
we observe that the cosmological term is a decreasing function of time and it 
approaches a small value as time progresses (i.e. the present epoch), which 
explains the small and positive value of A at present (Perlmutter et al. 
Riess et al. [T7|EH|; Knop et al. HI] Garnavich et al. [El; Schmidt et al. EI])- 
Figure 4 clearly shows this behaviour of A as decreasing function of time. 

The expressions for p{t), G{t) and A(t) for dust universe (p = 0, p > 0), 
Zel'dovich universe (ji = p), and radiative universe (p = |) can be obtained by 
putting 7 = 0, 1, i in Eqs. (j^SJl and respectively. 



G{t) 



A 



p{t) = B{ait^ + a2t + a3)"3(i+^\ 

[2aft2 + 2aia2t + al + k - 2aia3]{ait'^ + a2t + asy^+^t^ 
47r(l + 7)5 

^ _ (1 + 3-/){{2ait + a2f + fc} + ^a^ioLit^ + + as) 
(1 + 7)(aii2 +a2t + a3)2 



(22) 



(23) 



(24) 
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Figure 1: The plot of scale factor S{t) vs time with parameters ai = 0.1, 
02 = 0.01 and as — 1.0 
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Figure 2: The plot of energy density p{t) vs time with parameters ai = 0.1, 
a2 = 0.01, aa = 1.0, 7 = 0.0 and B = 1.0 
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Figure 3: The plot of gravitational term G{t) vs time with parameters k = 1.0, 
ai = 0.1, a2 = 0.01, as = 1.0, 7 = 0.0 and B = 1.0 



1.4= 

1.2- 




Figure 4: The plot of cosmological term A{t) vs time with parameter k = 1.0, 
ai = 0.1, a2 = 0.01, as = 1.0 and 7 = 0.0 
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6 Solution in the Sinusoidal Form 

Let L{R) = , where (3 is constant. 

In this case, on integrating, Eq. H13|) gives the exact solution 
S = M sin(/3i) + N cos{f3t) + (3i, 



(25) 



where M, N and f3i are constants. 

By using the same techniques as described in previous section 4, we find the 
expressions for p{t), G{t) and A{t): 



p{t) = B[M sin (3t + N cos pt + (3i 



-3(1+7) 



(26) 



G{t) = 



[f3^{{AP + N^)- f3i (M sin I3t + N cos /?<)} + k] [M sin pt + N cos pt + ^i](i+3t) 



A 



47r(l +7)B 
(1 + 37) [P^ ( M cos Pt-N sin ptf 



■k] 



( 1 + 7) [M sin pt + N cos pt + Pif 
2P^ {M sin pt + NCOS Pt) 



(27) 



(28) 



(1 + 7) [M sin pt + N cos pt + Pi] . 
Since, in this case, we have many alternatives for choosing values of M, TV, P 




Figure 5: The plot of scale factor S{t) vs time with parameter M = 1.0, N = 1.0, 
7 = 0.0, P = 0.10 and Pi = 2.0 
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Figure 7: The plot of gravitational term G{t) vs time with parameter k = 1.0, 
M = 1.0, N = 1.0, B = 10.0, 7 = 0.0, /? = 0.10 and /3i = 2.0 
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Figure 8: The plot of cosmological term A{t) vs time with parameter k — 1.0, 
M 1.0, N ^ 1.0, B = 10.0, 7 = 0.0, f3 = 0.10 and /?i = 2.0 

pi, it is just enough to look for suitable values of these parameters, such that 
the physical initial and boundary conditions are satisfied. We are trying to find 
feasible interpretation and situations relevant to this case. 

In Figure 5 we plot scale factor versus time. It is observed here that at 
early stages of the universe, the scale of the universe increases gently and then 
decreases sharply, and afterwards it will oscillate for ever. 

In Figure 6 we plot energy density p{t) versus time and in Figure 8 we note 
the run of cosmological constant A(t) versus time. From these figures we con- 
clude that at the early stages of the universe the matter is created as a result 
of loss of vacuum energy and at a particular epoch it has started to oscillate 
for ever due to sinusoidal property. In Figure 7 we plot the G{t) versus time. 
From it we conclude that in the beginning at the early stages of the universe 
G{t) increases gently with time and at a particular epoch it decreases sharply, 
and afterwards it will oscillate for ever. It is worth to mention here that in this 
case oscillation takes place in positive quadrant which is physically relevant. 

The expressions for p{t), G{t) and A{t) for dust universe {p — 0,p > 0), 
Zel'dovich universe {p = p), and radiative universe {p — ^) can be obtained by 
putting 7 = 0, 1, i in Eqs. (|2Sl, (ETjl and respectively. 
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7 Discussion and Conclusion 



In this paper we have described a new class of FRW cosmological models of 
the universe with a perfect fluid as the source of matter by applying a variable 
deceleration parameter. The "constants" G and A are allowed to depend on the 
cosmic time. Generally, the models are expanding, non-shearing and isotropic 
in nature. 

The cosmological terms in the models in Sections 4, 5 are a decreasing func- 
tion of time and this approaches a small value as time increases (i.e. present 
epoch). The value of the cosmological "terms" for these models are found to 
be small and positive which is supported by the results from recent supernovae 
la observations (Perlmutter et al. ^B]; Riess et al. Knop et al. j37) : 

Garnavich et al. |T2|; Schmidt et al. [2]). The cosmological term in Section 6 
also decreases while time increases to a specific instant. During this period as 
we can understand from Figure 8, we will have enough matter creation to force 
the universe to oscillate for ever due to sinusoidal property of A. This means 
we always have annihilation and creation of matter permanently. At this point 
one more sentence may be added to our discussion and i.e. as the graphs for 
A and p in this case the explosion of the universe at the early stages of its cre- 
ation has been only a consequence of matter creation. Thus, the implications of 
time varying A and G are important to study the early evolution of the universe. 
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